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PREFACE 


This manual describes the 32=bit run time library (RTL) 
mathematical functions. 


Chapter 1 introduces the RTL mathematical functions. Chapter 2 
gives the needed information for all subprograms. Chapter 3 
details elementary, min and max, and miscellaneous’ functions. 
The appendixes present information on error messages and the 
subprograms documented in the manual. Arithmetic tables are also 
included. 


The RO3 revision of this manual PEovaces more information on the 
exponentiation routines. 


For a survey of the Run Time Library and more information on the 
construction of programs in the Run Time Library, refer to the 
following Perkin-Elmer publication: 


@ FORTRAN VII Run Time Library (RTL) Introduction and Overview 
Reference Manual, 
Publication Number 29-578. | 


Additional information applicable to the Run Time Library is 
contained in the following Perkin-Elmer publications: 


@ FORTRAN VII Run Time Library (RTL) Language Extensions 
Reference Manual, 
Publication Number 29~580 


@e FORTRAN VII Run Time Library (RTL) Real Time Extensions 
Reference Manual, 
Publication Number 29-581 


e FORTRAN VII Reference Manual, 
Publication Number 48-017 


The following manuals provide detailed information on related 
software: 


@® Common Assembly Language (CAL) Programming 
Reference Manual, 
Publication Number 29-640 
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OS/32 Library Loader Reference Manual, 
Publication Number 48-020 


OS/32 Operator Reference Manual, 
Publication Number 29-574 


OS/32 Programmer Reference Manual, 
Publication Number 29-613 
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CHAPTER 1 
INTRODUCTION 


This manual describes the mathematical subprogram Perkin-Elmer 
provides in the FORTRAN VII Run Time Library. Subprograms 
described include the Intrinsic and Basic External functions 
required by ANSI Standard X3J3/77 and the arithmetic routines 
which the compiler calls. See the FORTRAN VII Run Time Library 
Language Extensions Manual, Publication Number 29-580, for the 
subprograms which perform bit and byte manipulation and logical 
operations. 


This manual provides the following information for each 
subprogram: 


- The purpose served by the subprogram. 

» The types of arguments and values. 

The range of admissible arguments. 

»- The accuracy (where relevant) of the result. 
- Error messages and possible interrupts. 

- A brief description of the algorithm. 


au & Wh 
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The subprograms can be called by assembly language code. The 
user must be aware of which arguments are passed to subprograms 
and which values are returned. See the FORTRAN VII Run Time 
Library Introduction and Overview, Publication Number 29-578. 


Chapter 2 contains information pertinent to all of the 
subprograms. Definitions are given to fix the terminology. 
Accuracy of the different subprograms, ranges of admissible 
values, and the generation of error messages are described. 


The Run Time Library provides descriptions of the function 
subprograms. These descriptions include external functions 
called by explicit statements in FORTRAN, intrinsic functions 
which can be called explicitly or implicitly (e.g., type 
conversion), and the involution routines called as a result of 
compiler-generated code. The programs are grouped as follows: 


Elementary functions 


Square root: SORT,DSQRT,CSQORT,CDSQORT 

Exponential: EXP,DEXP, CEXP, CDEXP 

Trigonometric: SIN,DSIN,CSIN,CDSIN,COS,DCOS,CCOS,CDCOS, 
TAN, DTAN 

Logarithm: ALOG,ALOG10,DLOG, DLOG10, CLOG, CDLOG 

Inverse Trigonometric: ACOS,ASIN, ATAN,ATAN2,CACOS, DASIN, DATAN, 
DATAN2 

Hyperbolic: SINH,COSH,TANH, DSINH, DCOSH, DTANH 


Ha-O025° FOO. ROG t= 


Soe ee eS 


Min and Max functions: 


AMAX0 AMAX]. AMINO MINO MAX] 
AMIN] DMAX 1] DMIN1 MAX0 MINI 


Type conversion subprograms: 


DBLE FLOAT IDINT DFLOAT 
IPIX INT INT 2 SNGL REAL 


Miscellaneous functions: 
Remainder 
AMOD DMOD MOD 
Manipulation of complex numbers 


AIMAG DIMAG CMPLX CON JG DCONJG 


Absolute value 

ABS CABS DABS DCABS IABS 
Transfer of sign 

DSIGN ISIGN SIGN 
Positive difference 

DIM IDIM DDIM 


Truncation and rounding 


AINT DINT ANINT DNINT NINT IDNINT 
Double~precision product of single~precision factors 


DPROD 


Involution (raising to powers): 


eI XXI -IXXR eI XXD eI XXC »~ITXXCD 
« RXXI «RKXR -RXXD «RXXC «RXXCD 
« DXXI »DXXR »-DXXD »DXXC -DXXCD 
- CXXI »~CXXR ~CXXD ~CXXC ~CXXCD 
- CDXXI »CDXXR ~CDXXD e«CDXXKC »~CDXCD 


Complex arithmetic: 


«CXMPY »~CDMPY ~CXDIV »~CDDIV 
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CHAPTER 2 
GENERAL CONSIDERATIONS 


2-1 TERMS RELATED TO ACCURACY 


Several methods describe the accuracy of a computation in 
floating— point arithmetic. In this manual, the following terms 
are used: 


Absolute Error: is the absolute value of the difference 
between the true value of a desired number and 
the value as computed. 


Relative Error: is the ratio between the absolute error and 
the true value. 


Generally, the relative error is the more useful of the two and 
better adapted to floating-point representation of numbers. To 
refer to a result as accurate "to five decimal places" is 
actually saying that the relative error is 1.0E-5 or better. 


A statement of the relative accuracy of a computation is based on 
the assumption that the numerical data, such as the argument of a 
function, is perfectly accurate. 


In practice, numerical data generally has errors. Therefore, it 
is important to know the manner in which these errors can grow in 
the course of a computation. This is discussed for each function 
or group of functions under the heading "Propagation of Errors." 


Extended discussions are availabie. for /- accuracy and 
error-propagation qualities of floating point programs. The 


reader can consult: 


Fike, C.T., Computer Evaluation of Mathematical Functions, 
Prentice-Hall, Inc., 1968 


Hart, J.F., et al, Computer Approximations, Wiley, 1968 
Sterbenz, P., Floating Point Computation, Prentice-Hall 

ae ACCURACY OF FUNCTION SUBPROGRAMS IN GENERAL 

The subprograms in the Run Time Library can be classified by 


their use of integer or floating-point arithmetic and by the 
amount of computation they involve. 


Floating-point computations: the function subprograms described 
under Elementary Functions generally return values with a 
relative accuracy better than l1.0E=-7 for Single-precision 


functions and better than 1.0E-16 for double-precision functions. 
There are exceptions to this rule, however, in the case of 
Logarithms and trigonometric functions. See Paragraphs 3.2.3 and 
3.2.4 for details. 
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The floating-point involution subprograms also have ae relative 
accuracy of 1.0E-7 for single precision or 1.0E-16 for double 
precision. The other computational subprograms are accurate to 
a few units in the last significant bit (i.e., a relative error 
of 1.0E-7 or better in single precision and 1.0E-16 in double 
precision). The value 1-0E-7 also applies to the subroutines 
mentioned under Complex Arithmetic, Paragraph 3.7. An exception 
in the case of AMOD or DMOD is described under Remaindering, 
Paragraph 3.5.1. 


Type conversion functions are accurate to one unit in the least 
Significant bit. Those subprograms using only integer arithmetic 
are absolutely accurate. The same is true of subprograms which 
essentially do no arithmetic (such as AMIN1 and SIGN). 


Ranges: The Basic External and Intrinsic function subprograms in 
the Run Time Library give meaningful results for all values of 
their arguments and are subject to the following range 
restrictions: 

Square root: The argument must be positive or zero. 
Exponential: The argument x must satisfy 


-~65*aloq(16) < x < 63*alog(16), i.e., 
-~180.2<x<174.6 


Logarithm: The argument must be positive. 


Trigonometric: The argument of SIN, COS, or TAN must) 
have an absolute value no greater than 
16**6 = 16,777,216. The argument of 
DSIN, DCcOS, or DTAN must have an 
absolute value no greater than 
16**14 = 72,057,594,037,927,9 36. 


Inverse trigonometric: The argument of ASIN, ACOS, DASIN, or 
DACOS must be between -1.0 and 1.0. 


Hyperbolic: The argument of SINH, COSH, DSINH, or 
DCOSH must have an absolute value no 
greater than 175.0. 


Complex functions: Restrictions on the arguments of 
complex functions are described under 
the individual functions. 


Type conversion: The argument of IDINT, INT, or MIFIX 
must be between -—(2**31-1) and 2**31-1 
(2**31 = 2,147,483,648). The argument 
of INT2, must be between -32,768 and 
32,767. 
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Miscellaneous functions: The second argument of AMOD, DMOD, or 
MOD, must be nonzero. The functions 
DSIGN, and SIGN require the second 
argument to be non-negative when the 
first argument is zero. 


Involution routines: The integer, real, double-precision, 
and complex routines reject a zero 
base with a negative exponent. The 
real and doublie-precision routines 
reject a negative base. 


2.3 GENERATION OF ERROR MESSAGES 


Even if arguments satisfy the restrictions detailed under Ranges, 
they can lead to arithmetic underflow or overflow. When these 
conditions do occur, or when it is known that they will occur, 
the Run Time Library function subprograms cause an error message 
to be printed. The same happens if inadmissible arguments are 
passed to the subprogram. Details on specific error messages are 
given for each subprogram under Description of Function 
Subprogramse The error messages are of the form: 


ERROR n IN m AT a: a 
where: n identifies the diagnostic message, 
m is the name of the subprogram, 


a is the hexadecimal address of the point in the 
user’s program from which the RTL is called, and 


dis a diagnostic message. 


See Appendix 1 for operation of the error message facility and 
suggested response to error messages. 


2.4 EXTERNAL AND INTRINSIC VERSIONS OF FUNCTIONS 


Except for the involution and complex arithmetic routines, every 
mathematical function in the RTL exists in two versions, as an 
intrinsic function or as an external function. The external 
functions are identified by their purely alphabetical names. The 
arguments are passed via an argument list. The intrinsic 
functions have names preceded by an @é. They also accept 
arguments via an argument list. Some intrinsic functions have 
argument passage in registers. Their names are preceded by a 
dot. Hence, certain functions in the RTL have 3 different 
versions, e.g. COS, @COS and .cos. 
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CHAPTER 3 


DESCRIPTION OF FUNCTION SUBPROGRAMS 


3.1 IWNTRODUCTION 


This section contains information on each function subprogram and 
it covers these topics: the purpose served by each function, the 
range of argument values which it accepts, considerations of 
accuracy and propagation of errors, an explanation of run-time 
error messages, subprograms called by a function, anda brief 
description of the algorithm is included for the more complicated 
mathematical functions. 


When a function is in both the intrinsic and the external 
versions, the external version calls the intrinsic version to 
perform the actual computation. Therefore, the references to the 
function subprograms use the names of. the intrinsic versions. 
Usually, the algorithms for single precision and double precision 
are Similar. To describe them concisely, the function name and 
its arguments written in lower-case letters are generic names. A 
reference for nearly all algorithms is the book by Hart, cited 
earlier in this manual. The largest number representable in 
floating-point notation is approximately 16**63 or 4.0E75. This 
number (whether in single or double precision) is referred to as 
"OMEGA". Small values of the argument are frequently a special 
case and dependent upon the level of precision. The notation 
"eps" denotes 16**(-3) in single precision and 16**(-7) in double 
precision: , 

EPS=Y *3E100000° in single precision 

EPS=Y~3A100000%, Y*00000000°% in double precision. 


The type Integer will always mean either Integer*2 or Integer*4. 


See E LEME NTA RY FUNCTIONS 


3.2.1 Square Root 


Argument Value 
Y=SQRT (xX) | Real Real 
DY=DSQRT (DX) Double precision Double precision 
CWW=CSQRT (CZZ) Complex*8 Complex*8 
CDWW=CDSQRT (CDZZ) Complex*16 Complex*16 
Purpose: To compute the square root of an argument 


of the appropriate type. 


Range of arguments: For SQRT and DSQRT, the argument must 
not be negative. For CSQRT and 
CDSQRT, the argument may be any 
complex number. 


Accuracy: SQRT and CSQRT are accurate to seven 
decimal places or better for all. 
arguments. DSQRT and CDS QRT are 
accurate to 16 decimal places or 
better for all arguments. 
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Propagation of error: 


Error messages: 


Subroutines called: 


Algorithms: 


The relative error in SQRT (X) is 
approximately half the relative error in 
xX. The only source of instability is the 
case where the absolute error in X is 
sizable compared to xX. 


If a negative x is received, then SQRT 
and DSORT print error messages "NEGATIVE 
ARGUMENT", The square root of the abso- 
lute value of x is returned to the called 
program. 


-CSORT calls .CABS and ~-SQORT 
-CDSORT calls .CDABS and .DSORT 


-DSQRT and .SQRT use the same basic 
algorithm. If x=0.0, then no processing 
is performed. Tf x<0.0, then an error 
message is printed and x is replaced by 
its absolute value. 


Let x=(16**m)*u, _ 1/16<u<l 
Let m=2*nt+ip, ip=0 or l 
u=v* (.25**ig), ig=0 or l, 1/4<v<l 


A rational approximation, wa, to sqrt(v) 
is found first (Hart, formulas 0291 or 
0293). The approximation is improved by 
Heron’s formula: 


wa is replaced by (v/watwa)*.5 


In double precision, this formula is 
applied twice. 


Let k=2*ip-ig 
Then sqrt (x)=(16**n)* (2**k)*wa 


The algorithm for CSQRT and CDSQRT 
makes use of SORT and DSORT, 
respectively. Let sqrt (xt+tiy)=atib. 
If x=y=0, no processing is performed. 


If x=0, then a=sart(abs(y)/2) 
b=(y/abs(y))a 


If y=0, x>0, then a=sqrt(x) 
b=0 

If y=0, x<0, then a=0 
b=sqrt(-x) 

Let r=cabs(xtiy) 


If x>0, then a=sqrt((xtr)/2) 
b=y/ (2a) 


If x<0, then b=(y/abs(y) )sqrt((-x+r)/2) 
a=y/ (2b) 
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3.2.2 Exponential Functions 


Y= EXP (X) 
DY=DEXP (DX) 
CWW=C EXP(CZZ) 
CDWW=CDEXP(CDZZ) 


Purpose: 


Range of arguments: 


Accuracy: 


Propagation of errors: 


Error conditions: 


Subroutines called: 
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Argument Value 
Real Real 
Double precision Double precision 
Complex*8 Complex*8 
Complex*16 Complex*16 

To compute exp (x) where x is a 


floating=point number of the appropriate 
type. 


For EXP and DEXP the argument must lie 
between -65*alog(16) and 63*alog(16). 
The interval from -180.0 to +174.0 is 
within these bounds. For CEXP and 
CDEXP, the real part must be between 
-180.0 and +174.0 and the imaginary 
part must have absolute value less 
than 16**6. 


EXP and CEXP normally give results 
accurate to seven decimal places or 
better. When the imaginary part of a 
complex argument is very large, this 


accuracy degrades; see remarks in 
propagation of errors under 
Trigonometric Functions, Paragraph 
3.2.3. DEXP and CDEXP give results 
accurate to 16 decimal places or 
better. 


The relative error of EXP(X) is equal to 
the absolute error of X. 


Tf the argument of EXP or DEXP is greater 
than the maximum acceptable value, then a 
message to that effect is put out and the 
number OMEGA is returned. If the 
argument of EXP or DEXP is less than the 
least acceptable value, then a different 
message is put out and 0.0 is returned. 
The same messages are put out if the real 
part of the argument of CEXP or CDEXP is 
out of the acceptable range. If the 
imaginary part of the argument of CEXP or 
CDEXP is unacceptable, another message is 
put out. CEXP returns the value (0,0) if 
the real part of the argument is less 
than the lower limit or the imaginary 
part is too large. If the imaginary part 
is acceptable, but the real part is above 
the upper limit, a large complex number 
is returned. 


-CEXP calls .EXP, .SIN, and .COS 


~CDEXP calls .DEXP, .DSIN, and .DC%S 


Algorithm: l. EXP and DEXP. 
It abs(x)<eps, then 1+x is returned. 
TE abs(x)>256 then it is treated as an 
errore Otherwise, let y=x*r, where 
r=2.0/alog(2.0) 
y=ntf£, O<f<l 
Tf n>504 then x is above the upper limit; 
1£f n<-520 then x is’ below the lower 
limit. Otherwise, exp (x) can be 
evaluated without overflow as follows: 
Let z=2**(g/2) 
z is found by a rational approximation 
(Hart, formulas 1080 or 1082). 
Let n=8*mtk, 0<k<7 
s=2**(0.5*k+0.25) 
Ss is found from a tabie 
Finally: 
exp (x)=(16**m) *s*z 
2. CEXP and CDEXP 
Let the two parts of the argument be 
(xr, xi) 
If -180.0<xr<174.0, let u=exp (xr) 
Tf xr<-180.0, let u=0 
TF xxr>174.0, let u=OMEGA 
If abs (xi)>16**6, let s=c=0. 
Otherwise, s=sin(xi), c=cos (xi) 
Then: 
cexp (xi, xr)=(u*c, u*s) 
3.2.3 Trigonometric Functions 
Argument Value 
Y =SIN(X) Real Real 
Y =COS(X) Real Real 
Y = TAN (X) Real Real 
DY =DSIN(DX) Double precision Double precision 
DY =DCOS( Dx) Double precision Double precision 
DY =DTAN (DX) Double precision Double precision 
Cww =CSIN(C2Z2Z) Complex*8 Complex*8 
CDWW=CDS IN (CDZZ) Complex*16 Complex*16 
CDWW=CDCOS (CDZZ) Complex*16 Complex*16 
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Purpose: 


Range of arguments: 


Accuracy: 


Propagation of errors: 


Error conditions: 


Subroutines called: 
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To compute the trigonometric sine, 
cosine, or tangent (assuming radian 
measure) of an argument of the 


appropriate type. 


For SIN, COS, and TAN, the argument 
must have an absolute value less than 
16**6. For DSIN, DCOS, and DTAN_ the 
argument must have an absolute value 
less than 16**4. For CSIN, CDSIN, 
CCOS, and CDCOS, the argument must 
have a real part less in absolute 
value than 16**6 and an imaginary part 
between -175 and +175. 


Normally, SIN, COS, TAN, CSIN, and 
CcOoS return values which are accurate 
to seven decimal places or better. 
DSIN, CDSIN, DCOS, and CDCOS return 
values accurate to 16 decimal places 
or better. 


The absolute error in SIN(X), COS(X), 
DSIN(X), and DCOS(X) is no greater on 


absolute value than the error in X. The 
absolute error in TAN(X) or DTAN(X) is 
the absolute error in xX times 


COS(X)**(-2). Thus when COS(X) is small, 
TAN(X) is large and its evaluation can be 


unstable. For the complex functions, 
when the imaginary part of the argument 
differs from 0.0, the relative error in 
the value is proportioned to the absolute 
error in the argument. 


If the argument of SIN, COS, or TAN is 
greater in absolute value than 16**6, 
then an error message is put out and 0.0 
is returned. Tf the argument of DSIN, 
DCOS, or DTAN is greater in absolute 
value than 16**14, then the same error 
message is put out and 0.0 is returned. 
Tf the argument of TAN or DTAN is such 
that function evaluation would lead to 
division by zero, then a different error 
message is put out and OMEGA is returned. 
If the real part of the argument of CSIN 
or CCOS has absolute value greater than 
16**6, a message is put out and (0.0, 
0.0) is returned. If the real part is 
acceptabie, but the imaginary part is 
greater in absolute value than 175.09, 
a message is put out and a large compiex 
number is returned. 


-CSIN and ~CCOS call OLN, -COSs, 
~-SINH, and .COSH 


-CDSIN and .~CDCOS call -DSIN, -DCOS, 
-DSINH, and .DCOSH 


Algorithms: 


1. 


Real and double-precision functions 


If abs(x)<eps, then 
sin(x)=x, cos(x)=1, tan(x)=x 


If abs(x) is too large (see "Error 
Conditions"), then 0.0 is returned. 


Otherwise let 


ya=x*2*c where c=2/7 
Set ya=k+u where -ls<u<l. 


Depending on k, the evaluation of the 
function of x reduces to one of these 
forms: 


pm*¥sin (u* 7/4) 
pm*cos (u* 7/4) 
tan (u* 7/4) 
cot (u* 7/4) 


where pm=l or -1l 


Polynomial or fractional approximations 
are used (Hart, formulae 3040, 3603, 
3820, 3823, 4242, and 4245) to evaluate 
the above 4 formulas. If u=0 and 
cot(u* 7/4) is required, then an error 
message is put out and OMEGA is returned. 


Complex functions 
Let the argument be cz2=(x,y).- 


If abs (x)<16**6, let 
rc=cos(x), rs=sin(x), 
Otherwise rc=rs=0 


If abs(y)<175.0, let 
hc=coshly), hs=sinh(y) ; 


Tf y>175.0, let 
hc=OMEGA, hs=OMEGA; 


hc=OMEGA, hs=-OMEGA 


Then: 


ccos(cz)=(re*hc, =-rs*hs) 
csin(cz)=(rs*he, re*hs) 
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3.2.4 Logarithmic Functions 


Y=ALOG (X) 
Y=ALOG10 (X) 
DY=DLOG (DX) 


DY=DLOG1LO (DX) 
CWW=CLOG (CZZ) 


Argument 


Real 

Real 

Double precision 
Double precision 
Complex*8 


Value 


Real 

Real 

Double precision 
Double precision 
Complex* 8 


CDWW=CDLOG(CD2ZZ) Complex*16 Complex*16 
Purpose: To compute the logarithm of a  floating- 
point number of the appropriate type. 
ALOG10 and DLOG10 compute logarithms to 


the base 10; the others, 
(natural logarithms). 


to the base e 


Range of arguments: ALOG, ALOG10, DLOG, 


positive (non-zero) 


and DLOGI1O require 
arguments. CLOG 


and CDLOG require a non-zero complex 
argument. 
Accuracy: ALOG, ALOG10O, and CLOG return values with 
an absolute error of 1.0E-7 or less. 
DLOG and DLOG10O return values with an 
absolute error of 1.0E-16 or less. 
Propagation of error: The absolute error of ALOG(X) equals’ the 


relative error of X. The absolute error 
of ALOG10(X) equals the relative error of 
X multiplied by ALOG10(e)=.434. 


Error conditions: When the argument of a real or 
double-precision function is 0, an 
error message is printed and -OMEGA is 
returned. When the argument is 
negative, a different error message is 
printed and the logarithm of the 
absolute value is returned. When the 
argument of CLOG or CDLOG is (0.0, 
0.0) an error message is printed and 
(-OMEGA 0.0) is returned. 

Subroutines called: -ALOG1 branches to .ALOG, and .DLOGI1 

. branches to .DLOG. 
-CLOG calls .~CABS, .ALOG1, and .ATAN2. 
-CDLOG calls -CDABS, -DLOGI, and 
-DATN2. 
Algorithms: 1. ALOG, ALOG10, DLOG, DLOG1O. 


After stage (A) in the algorithm, the 


computation for base 10 is identical to 
that for base e. Therefore, -ALOG1 
passes control to .ALOG, and .DLOGL 
passes control to .DLOG. 
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{(A) The error conditions are checked. 
The constant c is set to alog(2) in the 
programs .ALOG and .DLOG and to alogl0(2) 
in the programs .~ALOG] and .DLOGI1. 


(B) Let the argument be x 
x=(16**m)*£, 1/16<£f<l 
Let g=(2**j)*£, 0<3<3, .5<g<l 


LE g>sqrt(.5), let 
n=4*m=j 
h=(g-1)/(gt1); 


If g>sqrt (.5), let 
n=4¥*m-jel 
h=(g-.5)/(g+.5) 


The rational function of h (Hart’s 
formula 2701 or 2705) ais found; the 
result is the logarithm to the base 2 of 
either g or 2*g, depending on whether g 
is greater or less than sqrt(.5). tf 
this number is u, then v=ntu is the 
logarithm of x to the base 2s 
Multiplication by c gives the logarithm 
of x to the desired base. 


CLOG and CDLOG 


The error case (0.0, 0.0) is eliminated. 
If the argument (x,y) is not zero, let 


u=alog (cabs (x,y)) 
v=atan2 (y,x) 


Then (u,v) is returned 


36260 Inverse Trigonometric Functions 


Y=ACOS (X) 
Y=ASIN(X) 

Y=ATAN (X) 

Y=ATAN2 (Xl, X2) 
Y=DACOS (xX) 
Y=DASIN (X) 
Y=DATAN (X) 
Y2=DATAN2 (Xl, X2) 


ARGUMENT VALUE 

Real Real 

Real Real 

Real Real 

Real Real 

Double precision Double precision 
Double precision Double precision 
Double precision Double precision 
Double precision Double precision 
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Purpose: 


Programming note: 


Range of arguments: 


Accuracy: 


Propagation of error: 


Error conditions: 


Subroutines called: 


48-925 FOO P00 


To compute the inverse sine, cosine, or 
tangent of a number of the given type. 
The routines ATAN2 and DATAN2 compute the 
angle between an upward vertical line and 
the point with coordinates (x1l,x2). When 
x2>0, this is atan(x1/x2). In general, 
the value of ATAN2 is between -7 and 7. 


The definition of the value of ATAN2 is 
that which the ANSI standard requires. 
To get a value between 0 and 2* 7, one 
can take: 


gQMtatan2(-xl, -x2) when x1<0 


ATAN and DATAN accept any argument of the 
correct type. ATAN2 and DATAN2 accept 
any pair of arguments other than (0.0, 


0.0). ASIN, ACOS, DASIN, and DACOS 
accept any arguments between -1.0 and 
1.0 e 

The single-precision functions are 
accurate to seven places. The 


double-precision functions are accurate 
to 16 places. 


Both absolute and relative errors are 
decreased by ATAN and DATAN. When X is 
large, the absolute error in ATAN (X) is 
very small. For ATAN2 and DATAN2, it is 
convenient to define the relative error 
of a pair of arguments as: 


sqrt ((erl**2+er2**2)/(x1l**2+x2**2) ) 


where erl and er2 are the errors in xl 
and x2. Then, the absolute error in 
ATAN2 or DATAN2 is no greater than the 
relative error in the arguments. The 
functions ASIN, ACOS, DASIN, and DACOS 
increase the relative error; they become 
unstable as the argument nears 1.0 or 
=150.% 


If ATAN2 and DATAN2 receive (0.0, 0.0) as 
arguments, they put out an error message 
and return +7/4. If ASIN, DASIN, ACOS, 
and DACOS receive an argument greater 
than 1.0 in absolute value, they put out 


an error message and return the value 
0.0. 


-ASIN and .~.ACOS call .SQRT. e-DASIN and 
e-DACOS call .~DSQRT. 


Algorithms: 


1. 


ASIN, ACOS, DASIN, DACOS 


If abs(x)<eps, then 
asin(x)=x, acos(x)=7 /2-x 


If -.5<x<.5, then asin(x) is found by a 
rational approximation (Hart, 
formula 4693 or 4698); 
acos (x)=7 /2<asin (x). 


y=sqrt (0.5*x+0.5) 
u=asin (y) 

asin (x)=- 7 /2+2*u 
acos (x)= 7 /2=-2*u 


If e5<x<l, let 


y=sqrt (-0.5*x+0.5) 
u=asin(y) 

asin (x)= 7/2=-2*u 
acos (x)=2*u 


ATAN, DATAN 


If the argument is so large that the 
value of the function is +7/2 to the 
limit of precision involved, then that 
value is returned. Otherwise, the range 
is reduced by the following process. The 
argument belongs £o one of seven 
intervals, the dividing points being: 


-tan(75°), -tan( 45°), -~tan(15°), 
tan(15°), tan( 45°), tan(75°). The 
argument is shifted to the range 
~tan(15°)<w<tan(15°) by setting: 


w= (x—-a)/ (l+ax) 


where a=tan((30°)m), with -3<m<3 (in the 
extreme cases, the transformation is 
w=-1/x). Then atan(w) is evaluated; it 
is equal to w if abs(w)<eps, and 
otherwise w*P(w*w) where P is a 
polynomial. The coefficients of P are 
taken from Hart (No. 4940 for ATAN, No. 
4945 for DATAN). Then, 


atan (x)= atan(w)+m* (7/6). 
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The algorithms for ATAN2 and DATAN2 are 
like thase for ATAN and DATAN with the 
exception of the following: the case 
(0.0,0.0) generates an error message as 
mentioned above; in the case (xl, 0.0), 
the value +7/2 is returned if x1>0.0 and 
- W/2 if x1<0; in the case (0.0,x2), the 
value 0.0 ais returned if x2>0 and Wif 
x2<0. 


The exponents of the floating=point 
numbers xl and, x2 are checked; and, if 
they differ by more than 6, the smaller 
is treated as 0.0 according to the 
previous paragraph. The exception is the 
case when xl=0.0 and x2>0; then xl/x2 is 
returned. 


In the general case, the argument pairs 
are partitioned into sectors depending on 
the signs of the arguments and the size 


of their ratio. The evaluation is 
reduced toa finding atan (w) where 
w<tan(15°). The value of atan(w) is 


computed exactly as in ATAN or DATAN. 


3.2.6 Hyperbolic Functions 
- ARGUMENT VALUE 
Y=SINH (Xx) Real Real 
Y=COSH (Xx) Real Real 
Y= TANH (X) Real Real 
Y=DSINH(X) Double precision © Double precision 
Y=DCOSH (X) Double precision Double precision 
Y=DTANH(X) Double precision Double precision 


Range of arguments: 


Accuracy: 


Propagation of error: 
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SINH, COSH, DSINH, and DCOSH accept 
arguments between -175.0 and +175.0. 
TANH and DTANH accept all values. 


The single-precision functions are 
accurate to seven decimal places; the 
double-precision functions are accurate 
to 16 decimal places. 


SINH multiplies the absolute error in xX 
by COSH(X) and COSH multiplies it by 
SINH(X). When xX is greater than 1 in 
absolute value, the relative error in 
SINH(X) or COSH(X) is close to . the 
absolute error in XX. TANH decreases both 
relative and absolute errors. 


Error conditions: 


Subroutines called; 


Algorithms: dng 


S=b2 


If the argument is greater than 175.0, 
then SINH and cCOSH put out an error 
message and return OMEGA. if the 
argument is less than “175.0, then SINH 
and COSH put out an error message and 
return, respectively, sOMEGA and OMEGA. 


-SINH and ~-COSH call .HYEXP and - TANH 
calls .EXP. 


~-DSINH and .DCOSH call .HYDEX, and .DTANH 


-Calls .DEXP. 


SINH and DSINH 
If abs(x)<eps then sinh(x) =x 


If abs(x)>175.0 then an error condition 
exists as described above. 


Lf -175.0<x<-.5, let 


xa=-x-alog(2) 
=-exp (xa) 


If u<-16**3 in single precision or 
u<-16**7 in double precision, 


then’ 
sinh (x) =u, 
otherwise 
sinh(x)=u-.25/u 


Tf -.5<x<.5, a rational approximation is 
used (Hart, formula 2023 or 2028). 


If .5<x<175.0,let 


“xa=x-alog (2) 
u=exp (xa) 


if u>lo**? 17 single precisicn oF 
u>l6**7 in double precision, 


then 
sinh(x) =u, 
otherwise 


sinh(x)=u-.25/u 
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2. COSH and DCOSH. 
If x<0, replace x by -x 
Lf x>175.0, an error condition exists as 


described above. If x<eps, 


Then 
cosh(x)=1.0 


Otherwise let 


xa=x~alog(2) 
u=exp (xa) 


If u>lo**3 in single precision or u>16**7 
in double precision, 


Then 
cosh(x)=u 
Otherwise 
cosh(x)=ut+.25u 
3. TANH and DTANH 


Let top=8.5 in single precision and 
21.0 in double precision. 


If x>top, tanh(x)=1 
If x<-top, tanh(x)=-1 
If .5*alog(2)<x<top, let 


u=exp (2*x) 
tanh (x)= (u-1)/ (utl1) 


If abs(x)<eps, then tanh(x)=x 


If eps<abs(x)<.5*alog(2), a rational 
approximation is used. 


If -top<x<-.5*alog(2), let 


u=exp (-2*x) 
tanh (x)=(l-u)/ (itu) 
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Purpose: 


MIN AND MAX FUNCTIONS: 


VALUE 


ARGUMENT 
I=MaAxX0O (Il,...,I1n) Integer Integer*4 
I=MINO (11,...,I1n) Integer Integer *4 
X=AMAXO (1I1,...,In) Integer Real 
X=AMINO (I1,...,In) Integer Real 
I=MAX1 (Xl,...,Xn) Real Integer*4 
I=MIN1 (Xl,...,Xn) Real Integer*4 
X=AMAX1 (Xl,...,Xn) Real Real 
X=AMINI1 (X1l,...,Xn) Real Real 
DX=DMAX1 (DXl,...,DXn) Double Double 
DX=DMINI1 (DX1l,...,DXn) Double Double 
where 2<n<256 
To find the algebraic minimum or maximum 
of a collection of arguments of a given 
type and (in some cases) to convert toa 
different type. Types of arguments) and 
values are shown in Table 3-l. 
TABLE 3-1 TYPES OF MIN AND MAX ARGUMENTS 


INTEGE 


DOUBLE 


R*4 


PRECISION 


ARGUMENT 


DOUBLE 


INTEGER REAL 


PRECISION 


MAX 0 
MINO 


AMAX 0 
AMINO 
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Range of arguments: 


Accuracy: 


Propagation of error: 


Error messages: 
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Unrestricted, except for MAX] and MINI, > 
which convert real numbers to integers. 
MAX] and MINIL require the chosen value to 
be between -2,147,483,648 and 
+2,147,483,647. 


Absolute. 


Any error in the argument chosen as 
minimum or maximum is passed through 
unchanged. 


If MAX1 or MINI, try to convert a number 
too large for the appropriate type to 
fixed point, they print “OVERFLOW and 
return the largest integer of the 
appropriate sign and type. 


3.4 TYPE CONVERSIONS 


Purpose: 


Range of arguments: 


Accuracy: 


Propagation of error: 


Error messages: 


Remarks: 


To convert a number from one type to 
another. The types of arguments and 
values are shown in Table 3-2. 


The functions with INTEGER*2 value 


require arguments between ~-32,768 and 
32,767. The functions with INTEGER*4 
value require arguments between 


-2,147,483,648 and +2,147,483,647. 


SNGL truncates the least significant 32 
bits of the fraction in a 
double~precision number. The functions 
which convert floating point to integer 
all truncate those bits which represent 
values less than one. For large 
arguments, FLOAT also truncates some 
bits; the result is accurate to the least 
significant bit in a floating point 
number. 


A small error in the argument can be 
entirely lost or it can result ina large 
error in the value, depending on whether 
or not it moves the argument from below 
an integer to above that integer. If the 
true value of x is 2.93, then INT(X)=2. 
If xX ais computed as 2.85 (an error of 
-08), then INT(X) is correctly computed 
as 2. But if X is computed as 3.01 (also 
an error of .08), then INT(X) is computed 
as 3. 


If the functions with integer values 
receive an argument which is too large, 
they print "OVERFLOW" and return the 
largest integer of the appropriate sign 
and type. 


The name INT2, INT, REAL, DBLE, and CMPLX 
are generic names for FORTRAN VII O and 
FORTRAN VII _ OD. This means that INT2, 
INT, REAL, DBLE, and CMPLX can take any 
of the following types of arguments: 


- Integer*2 

- Integer*4 

- Real*4 

- Double Precision (Real*8) 
- Complex*8 

- Complex*16 
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Table 3-2 lists their names, types of arguments, and return 
values. 
TABLE 3-2 TYPE CONVERSION ROUTINE NAMES, TYPES OF 
ARGUMENTS AND VALUES 
TYPE OF INPUT 
Specific | FORTRAN ARGUMENT 2 TYPE 
Generic RTL User Optional OF 
Name Name Name Arg l | Arg 2 FUNCTION 
Integer Integer*2 
INT2 Real Integer*2 
Double Integer*2 
Complex*8 Integer*2 
Complex*16 Integer*2 | 
INT Integer Integer 
r 
eee Real Tateges 
IFIX 
LDINT Double | Integer 
Complex*8 Integer 
Complex*16 
REAL Integer 
@F LOAT Integer 
@@RREA Real 
@SNGL Double 
@@CREA Complex*8 | 
@@CDRE Complex*16 
: @DFLOAT Integer 
ale @DBLE a Real 
@@CDBL Double 
@@DDBL Complex*8 
@DREAL Complex*16 . 
CMPLX Integer In teger Complex*8 
CMPLX Real Real Complex *8 
Double Double |} Complex*8 
Complex*8 Complex*8 Complex *8 
Complex*16 oe A _| Complex*16_ 
| pempLx @@IcpD Integer Integer | Complex*16 
@@RCD | Real*4 Real*4 Complex*16 
@DCMPL DCMP LX Real*8 Real*8 Complex*16 
@@CxCD Complex*8 Complex*l16 
| @@CDCD Complex*16 Complex*l 6 
CMPLX or DCMPLX takes 2 arguments of 
the same type, if the type is not 
complex. CMPLX or DCMPLX will accept 
only 1 argument of the complex type, 
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in which case CMPLX just returns) that 


argument. 


A complex number CZ is stored as two real 
numbers (21, 22) signifying (Z1) + (22) * 
SORT(-1). The statement CZ COMPLX (Al, 
A2) for Al, A2 of the same type which is 
not complex, results in making Al, A2 
into 2 real numbers 21, 22 (by floating 
or truncation) and setting CZ C20 22). 


Soe a 


3.5 MISCELLANEOUS FUNCTIONS 


3.5.1 Remaindering 


Y=AMOD (Xl, X2) 
DY =DMOD (DX1, DX2) 
K=MOD (11, 12) 


Purpose: 


Range of arguments: 


Accuracy: 


Propagation of error: 


krror messages: 


ARGUMENT VALUE 

Real. Real 

Double precision Double precision 
Integer Integer *4 


To find the remainder’ as a result of 


dividing the first argument by the 
second. 
The second argument must not be zero. In 


the functions AMOD and DMOD, if the 
second argument is very small compared to 
the first, then a floating point overflow 
can occur. 


MOD iS’ precise. AMOD is accurate to 
about 1.0E-7 times the value of X2 and 
DMOD to about 1.0E-16 times the value of 
DX2. These accuracies can be degraded 
when X1/X2 becomes large. When X1/X2 is 
greater than 16**6, the returned result 
is meaningless. 


Absolute error in the first argument of. 
AMOD or DMOD is passed through unchanged, 
unless it changes the integer part of 
X1/X2. In that case, it leads to an 
error of the same size as X2. Absolute 
error in the second argument is 
multiplied by the integer part of Xl/X2; 
if it also causes a shift in that number, 
then an error the _ size of X2 is 
introduced. 


Examples: 


Let the true value of x be 10.0 and the 
true value of y, 91.0. Then amod (y,x)= 
1.0. If y is computed as 92.1, then 
amod(y,x)=2.1; the error is passed along. 
If y is computed as 89.9 (same error in 
opposite direction), then amod(y, x)=9.9; 
a large error is introduced. If y is 
computed as 91.0 and x is computed as 
10.1, then amod (y,x)=0.1; the error‘in x 
has been multiplied by y/x. If x is 
computed as 10.2, then amod (y,x)=9.4; a 
large error is introduced. 


If the second argument is 0, then an 
error message is printed and the value of 


the first argument is returned. 
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3.5.2 Manipulation of Complex Numbers 


_ARGUMENT VALUE 
Y=AIMAG(C2Z) : Complex*8 Real 
R= CABS (CZ) Complex*8 Real 
CW=CONJG(CZ) Complex*8 Complex* 8 
CDWW=DIMAG(CDZZ) Complex*16 Complex*16 
CDWW=CDABS (CDZZ) Complex*16 Complex*16 
CDWW=DCONJG(CDZZ) Complex*16 Complex*16 
Purpose: The complex number CZ is stored as two 
real numbers (Z1, 22) % signifying 


(Z1)+(Z2)SQRT(-1). AIMAG and DIMAG 
return 22. CABS and CDABS return 
SQRT((Z1)*(Z1)+(22)*(Z2)). CONJG and 
DCONJG return (21, -22) as its value. 


Range of arguments: Unrestricted. 


Accuracy: Not applicable except to CABS, which 
is accurate to seven decimal places; 
and CDABS to sixteen decimal places. 


Error messages: CABS and CDABS can cause arithmtic 
overflow. 


Subroutines called: CABS calls .SQRT 


CDABS calls .DSQRT 


3 ose Absolute Values 


ARGUMENT _ VALUER 
Y=ABS (X) Real Real 
DY=DABS (DX) Double precision Double precision 
T=IABS (J) Integer Integer *4 
Purpose: To return the absolute value of an 


appropriate argument. 


Range of arguments: Unrestricted, except LABS return 
+2,147,483,647 for full scale negative 
number -2,147, 483,648. 


Accuracy: Not applicable. 


Propagation of error: Not applicable. 


Error messages: Not applicable. 


Subroutines called: Wone.e 
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3.5.4 Transfer of Sign 


K=ISIGN(I,J) 
Z=SIGN(X,Y) 
DZ=DSIGN(DX,DY) 


Purpose: 
Range of arguments: 
Accuracy: 


Propagation of error: 


Remark: 


Error messages: 


Subroutines called: 


ARGUMENTS VALUE 
Integer Integer *4 
Real Real 


Double precision Double precision 


The second argument sign is applied to 


the absolute value of the first. 


the first 


The only restriction is 
if the second 


argument must not be 0 
argument is negative. 


Does not apply. 
Does not apply. 


zero is a positive number for these 


routines. 


If the first argument is 0 and the second 
argument is negative, then an error 
message is printed and 0 is returned. 


None. 


3.5.5 Positive Difference 
ARGUMENT VALUE 
K=IDIM(I,J) Integer Integer*4 
2=DIM(X,Y) Real Real 
DZ=DDIM(DX,DY) Double precision _ Double precision 
Purpose: To give the value 0 or the value of the 
first argument minus the second, 
whichever is greater. 
Range of arguments: Unrestricted. 
Accuracy: To one bit in the least significant 
place. 
Propagation of error: Not applicable, 
Error messages: Not applicable. 
Subroutines called: None, 
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Y=AINT(X) 
Y=ANINT(X) 
DY=DINT (Dx) 
DY=DNINT (DX) 
J=NINT(X) 

J= IDNINT (DX) 


Truncation and Rounding 


Argument 


Real 
Real 
Double precision 
Double precision 
Real 
Double precision 


Value 


Real 

Real 

Double precision 
Double precision 
Integer 

integer 


Function 


Truncation 
Rounding 
Truncation 
Rounding 
Rounding 
Rounding 


To obtain integer values with truncation, use the type conversion 
routines IDINT,.INT, or INT2. 


“The Truncation functions remove the 
fractional part of a number; the result 
is equal to or smaller in absolute value. 
The Rounding functions return the nearest 
integer value. A fractional part of .5 


Purpose: 


is rounded up in absolute value; thus, 
Range: Unrestricted for AINT, ANINT, DINT, and 


DNINT. 


For NINT and IDNINT, the value must be no 
greater than 2**31-1 in absolute value. 


Accuracy: Not applicable. 


Propagation of error: Not applicable. 


Error messages: An overflow condition can be incurred by 


NINT or IDNINT. 


3.5.7  Double-Precision Product of Single-Precision Factors 
Argument Value 
DY=DPROD(X1,X2) Real Double-precision 


Purpose: The argument’s product is computed as a 


double-precision number. 


Range of Arguments: That of floating-point multiplication. 


Accuracy: No error is introduced. 


Propagation of error: the 


the 


The relative error in the product is 
sum of the relative errors in 
factors. 
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Error messages: 


Subroutines called: 


3.6 INVOLUTION (RAISING 


Program 


- TXXI 
» TXXR 
~-IXXD 
e IXXC 
» TXXCD 
e RXXI 
~-RXXR 
-RXXD 
« RXXC 
« RXXCD 
- DXXT 
- DXXR 
- DXXD 
« DXXC 
- DXXCD 
« CX XI 
« CXXR 
~CXXD 
« CX XC 
» CXXCD 
« CDXXI 
«CDXXR 
« CDXXD 
« CDXXC 
- CDXCD 


Purpose: 


Range of arguments: 


Multiplicative 


overflow or 


underflow 


results in error messages. 


Integer 

Integer 

Integer 

Integer 

Integer 

Real 

Real 

Real 

Real 

Real 

Double Precision 
Double Precision 
Double Precision 
Double Precision 
Double Precision 
Complex 

Complex 

Complex 

Complex 

Complex 

Complex *16 
Complex* 16 
Complex *16 
Complex*16 
Complex*]1 6 


To 
a**b, 


None. 


evaluate 


TO POWERS) 


Ex ponent 


Integer 

Real 

Double Precision 
Complex 
Complex*16 
Integer 

Real 

Double Precision 
Complex 
Complex*16 
Integer 

Real 

Double Precision 
Complex 
Complex*16 
Integer 

Real 

Double Precision 
Complex 
Complex*16 
Integer 

Real 

Double Precision 
Complex 
Complex*16 


(for 


defined as the 


determined by 
arithmetic 


bases and exponents. 


of 


exponentiation 
mode data 


types is 


expressions of the 

complex involutions a**b is 
principal 

EXP (b*LOG (a) )). 

types are 


Value 


Integer*4 


Real 

Double Precision 
Complex 
Complex*16 

Real 

Real 

Double Precision 
Complex 
Complex*16 
Double Precision 
Double Precision 
Double Precision 
Complex*]6 
Complex*16 
Complex 

Complex 

Complex *16 
Complex 
Complex*16 
Complex*16 
Complex*16 
Complex*16 
Complex*16 
Complex*1 6 


form 
value 


All 
permissable as 


The interpretation 
operations 
described in the 


on mixed 


FORTRAN VII Reference Manual. 


These 
acceptable: 
exponent; 


following 
zero 


exponent. 


conditions are 
base 
negative base with nonwinteger 


not 


with negative 
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Accuracy: 


Propagation of error: 


Interfaces: 


Error conditions: 
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~-IXXI 
- RXXR, 
decimal 
~CDXXI 
places. 


The relative 
relative error of the base multiplied by 
and the absolute error of 
natural 


the exponent 
the exponent multiplied by 


and 


introduces 
«CXXI 


places. 


no 


-DXXI, 


are accurate to 


error 


logarithm of the base. 


Arguments and values. are 
loading them into registers as follows: 


Program 


eI XXI 
-ITXXR 
~IXXD 
e TXXC 
-IXXCD 
« RXXI 
«RXXR 
e RXXD 
« RXXC 
e RXXCD 
» DXXI 
« DXXR 
- DXXD 
« DXXC 
- DXXCD 
« CX XI 
« CXXR 
« CXXD 
« CXXC 
« CX¥XCD 
- CDXXI 
- CDXXR 
« CDXXD 
« CDXXC 
« CDXCD 


F14 
F14 


F14 


D14 

D14 

D14 

D14 

D14 
(F12,F14) 
(F12,F14) 
(F12,F14) 
(F12,F14) 
(F12,F14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 


Exponent 


Gl2 

F14 

D14 
(F12,F14) 
(D12,D14) 

G13 

F12 

D14 
(F10,F12) 
(D12,D14) 
G13 

F14 

D12 
(F12,F14) 
G13 

F10 

D14 
(F8,F10) 
(D12,D14) 
G13 

F14 

D10 
(F12,F14) 
(D8,D10) 


Unacceptable combinations of 
are detected and error messages are put 
out. For a zero 


exponent, 
returned. 


base and 


error. 
are accurate to six 
e-DXXD, 


includes: 


the 


passed 


eRXXI, 


decimal 


Value 


D14 
(F12,F14) 
(D12,D14) 

F14 

F14 

D14 
(F12,F14) 
(D12,D14) 
D14 

D14 

D14 
(D12,D14) 
(D12,D14) 
(F12,F14) 
(F12,F14) 
(D12,D14) 
(F12,F14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 
(D12,D14) 


arguments 


negative 


the largest positive number is 


For 


floating~point 


returned. 


and -DXXD can 
occurrences of arithmetic underflow or 


overflow; 
occurrences 
other than 


eI XXI 


negative 


exponent, 
The routines 


detects 
internally. 
eIXXI can suffer 


The 
exponential 


base and 


0.0 is 


«IT XXI, .~RXXR, 
anticipate 


certain 


all such 


routines 


and 


the 


by 
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Algorithm: 


Remarks: 


underflows or overflows that cause 
arithmetic fault interrupts. In these 
cases, the user sees error messages from 
both the operating system and the Run 
Time Library. 


If either a fixed-point overflow or a 
floating point exponential overflow is 
detected, the value returned is the 
largest number of the appropriate sign. 
If an exponential underflow is detected, 
the value returned is 0.0. 


If the exponent is of integer type, the 
subroutine repeatedly squares the base 
and multiplies together those iterated 
squares that correspond to “1” bits in 
the binary expression of the exponent. 


The routines which receive real 
exponents of a single or double 
precision use the formula 


x**y=exp(y*alog(x)). If the exponent is 
0, then the result 1 or 1.0 is 
invariably returned, whatever the value 
of the base. 


These routines are called by the FORTRAN 
compiler to evaluate expressions of the 
form a**b. User can access’ these 
routines via assembly code. 
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ERROR MESSAGES 


If a mathematical RTL program detects an error, it calls a 
routine ~.ERR to print out an error message in this form: 


ERROR n IN m AT a; d 


where: nis the error message number, 


m is the name of the RTL function, 


a is the hexadecimal return address in the 
user”’s program, and 


ad is either a number or a diagnostic message. 


The messages are listed by number in the error messages list in 
this appendix. The following comments are intended to help 
interpret the messages and correct their causes. 


Messages Pertaining to the Argument List 


When a mathematical RTL program is called by giving it an 
argument list with argument checking turned on, ae series of 
messages can be put out as a result of errors in the list. 


Please refer to FORTRAN VII Run Time Library Introduction and 
Overview, Publication Number 29-578 for argument checking. 


The error messages may be one of the following: 


1. ILLEGAL CALL (72). 
A function has been called as a subroutine, e.g. 


CALL SQRT (cece) 


The solution is to revise the text of the source 
program or, if the user means to call an external 
subroutine of that name, to ensure that it is brought 


in at linkage time. 


2e INCORRECT TYPE OF VALUE OF FUNCTION (73). 
A function supplied in the RTL has been called as an 
external function and the user program has specified a 
value type that is incorrect. Example: 


IMPLICIT INTEGER (A~-H) 
EXTERNAL COS 
A=COS (X) 


Check implicit and explicit types or include a 
user=written program of the appropriate name. 
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3. ARGUMENT OF INCORRECT CLASS (75). 
An argument has been classified aS an array or 
external function and not a number. Check for the 
presence or absence of a DIMENSION statement, a 
misnamed argument, or an incorrect EXTERNAL statement. 


4. AN ARGUMENT OF INCORRECT TYPE (45). 
The argument supplied is of the wrong type for this 
function. Check implicit and explicit type statements 
and types of expressions passed as arguments. 


5. INCORRECT NUMBER OF ARGUMENTS (76). 
Check argument lists in source program and for the 
presence or absence of commas. 


Arguments Outside Acceptable Ranges 


These messages can occur whether or not a program is called with 
arguments in registers or by giving it a list. 


1. ARGUMENT TOO LARGE (49) 


REAL PART TOO LARGE (52) 
IMAGINARY PART TOO LARGE (55) 


The absolute value of the number is too large if the number is 


positive or negative. Check the numerical stability of the 
computation being performed. Check if variables have been 
initialized. Consider building errorechecking logic into the 
program. 


2. ARGUMENT OUT OF RANGE: POSITIVE (50) 
ARGUMENT OUT OF RANGE: NEGATIVE (51) 
REAL PART OUT OF RANGE: POSITIVE (53) 
REAL PART OUT OF RANGE: NEGATIVE (54) 
IMAGINARY PART OUT OF RANGE: POSITIVE (56) 
IMAGINARY PART OUT OF RANGE: NEGATIVE (57) 


The number is greater than its maximum permissible value or less 
than its minimum permissible value. Remedies are as in #1. 


3. OVERFLOW ON CONVERSION (58) 
The absolute value of a floating-point number was too large to be 


converted to an integer. Besides the messages listed in 1 above, 
consider if this conversion is really appropriate. 
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4. ZERO ARGUMENT (46) 
NEGATIVE ARGUMENT (47) 
ARGUMENTS (0,0) (48) 
ZERO DIVISOR (62) 
ARGUMENTS (ZERO, NEGATIVE) (63) 


These cases must be averted by building appropriate logic into 
the source code. Numerical instability should be closely 
checked; the difference of two numbers may be quite inaccurate. 


5. NEGATIVE EXPONENT (64) 
NEGATIVE BASE (74) 
ZERO BASE, NEGATIVE EXPONENT (65) 


These are conditions that can arise in computing powers. Check 
the validity of the computations which determined the base and 
exponent; consider building in logic to avoid these erroneous 
conditions. 


6. EXPONENTIAL OVERFLOW (60) 
EXPONENTIAL UNDERFLOW (61) 


These are conditions that can arise when the computation of a 
number raised to a power has ted to arithmetic faults. To 
safeguard against these conditions, build logic into the program 
to handle extreme values of either base or exponent. 
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RTL PROGRAM ERROR MESSAGES 


ERROR MESSAGES THAT MAY BE INVOKED FROM MATHEMATICAL RTL 
PROGRAMS. 


NUMBER DIAGNOSTIC MESSAGE 
45 AN ARGUMENT OF INCORRECT TYPE 
46 ZERO ARGUMENT 
47 NEGATIVE ARGUMENT 
48 ARGUMENTS (0,0) 
49 ARGUMENT TOO LARGE 
50 ARGUMENT OUT OF RANGE : POSITIVE 
51 ARGUMENT OUT OF RANGE : NEGATIVE 
52 REAL PART TOO LARGE 
53 REAL PART OUT OF RANGE : POSITIVE 
54 REAL PART OUT OF RANGE : NEGATIVE 
55 IMAGINARY PART TOO LARGE 
56 IMAGINARY PART OUT OF RANGE : POSITIVE 
57 IMAGINARY PART OUT OR RANGE : NEGATIVE 
58 OVERFLOW ON CONVERSION 
59 VALUE OF SIZE IS ILLEGAL 
60 EXPONENTIAL OVERFLOW 
61 EXPONENTIAL UNDERFLOW 
62 ZERO DIVISOR 
63 ARGUMENTS (ZERO, NEGATIVE) 
64 NEGATIVE EXPONENT 
65 ZERO BASE, NEGATIVE EXPONENT 


Al-4 48+~025 FOO ROO 


NUMBER 


APPENDIX 1 (Continued) 


RTL PROGRAM ERROR MESSAGES (Continued) 


66 


67 


68 


69 


70 


rank 


72 


13 


74 


75 


76 


DIAGNOSTIC MESSAGE 
NEGATIVE INDEX 
SIZE TOO LARGE 
SIZE NOT POSITIVE 
ILLEGAL FUNCTION CODE 
ILLEGAL TRAP CODE 
ILLEGAL LEVEL OF TRAPPING 
ILLEGAL CALL 
INCORRECT TYPE OF VALUE OF FUNCTION 
NEGATIVE BASE 
ARGUMENT OF INCORRECT CLASS 


INCORRECT NUMBER OF ARGUMENTS 
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APPENDIX 2 
ALPHABETICAL LIST OF SUBPROGRAMS 
DOCUMENTED IN THIS MANUAL 


Math Paragraph _Math Paragraph 
Function Reference Function Reference 
@ABS 3.5.3 @DFLOAT 3.4 
@ACOS 3.2.5 @DIM 3.5 6 
@AIMAG 3.5.2 @DIMAG 3.5.2 
@AINT 3.5.6 @DINT 3.5.6 
@ALOG 8.6254 @DLOG 22224 
@ALOGL 34 Or @DLOGL 330.4 
@AMAX0 3.3 @DMAX1 3.3 
@AMAXL 3.3 @DMIN1 3.3 
@AMINO 3.3 @DMOD 3.5.1 
@AMIN1I 3.3 @DNINT 34546 
@ANINT iris ae @DPROD S467 
@ASIN 3.2.5 @DREAL ,3-4 
@A TAN 320.5 @DSIGN 3.5.4 
@ATAN2 3.2.4 @DSIN 3.2.3 
@DSINH 3.2.6 
@DSQRT Bi 2 chp Ss 
@CABS 352-54 @DTAN a223 
@ccos 359) 3 @DTANH 3.2.6 
@CDABS 23522 
@cpcos 3.253 @EXP 3:32 
@CDEXP 3.2.2 @FLOAT | 3.4 
@CDLOG 3.2.4 @IABS 3.523 
@CDSIN 262.3 @IDIM S55 
@CDSQR $3.4 @IDINT 3.4 
@CEXP 320-2 @LDNIN 3.5.6 
@CLOG oe @IFIX 3.4 
@CMPLX 3.4 @INT 3.4 
@CONTG 3.5.2 @INT2 3.4 
e@cos 349.5 @ISIGN 3.5.4 
@COSH 32:3 
@CSIN a7 955 @MAX0 3.3 
@CSQRT cee ee @MAX1 33.3 
@MINO 3.3 
@MINL 2.3 
@DABS 3.5.3 @MOD 3.5.1 
@DACOS 3.2.5 @N INT 3.5.6 
@DASIN 3.2.5 
@DATAN 3.2.5 @REAL 354 
@DATN2 3.2.5 @SIGN 3.5.4 
@DBLE 3.4 @SIN ay Dos S 
@DCMPL 34:4 @SINH Seah Bis 
@DCONJG 34562 @SNGL 3.4 
@Dcos 3.2.3 @SQRT 3.2 
@DCOSH 3. D6 
@DDIM 3.5.5 @TAN a 2a 
@DEXP Sno 6 @TANH 3.2.6 


43-925 FOO ROO 


e 
bd 
. 
~ 
W 
. 
nN 
e 
fon) 


wow ee Be ow 


woe we ee we 


Math Paragraph Math Paragraph 


Function Re ference Function Re ference 
@@CCX 3.4 - DACOS 3.2.5 — 
@@C DBL 3.4 -DASIN 3.2.5 
@@CDCD 3.4 - DATAN 3.2.5 
@@CDCX 3.4 - DATN 2 3.2.5 
@e@cpDI 3.4 - DCOS 3.2.3 
@@CDI2 3.4 - DCOSH 3.2.6 
@@C DRE 3.4 - DEXP 3.2.6 
@@CI2 3.4 -DINT 3.5.6 
@@CINT 3.4 - DLOG 3.2.4 
@@CREA 3.4 -DLOG1 3.2.4 
@@CXCD 3.4 -DNINT 3.5.6 
@@ICD 3.4 - DSIN 3.2.3 
@@RCD 3.4 -DSINH 3.2.6 

- DSQRT 3.2.1, 3.2.5 
-ACOS 3.2.5 - DTAN 3.2.3 
-AINT 3.5.6 - DTANH 3.2.6 
- ALOG 3.2.4 »DXXC 3.6 
- ALOG1 3.2.4 » DXXCD 3.6 
-ANINT 3.5.6 -«DXXD 3.6 
-ASIN 3.2.5 -DXXI 3.6 
-ATAN 3.2.5 -DXXR 3.6 
- ATAN2 3.2.4 

- EXP 3.2.2, 3.2.6 
-CABS 3.2.4 
- CCOS 3.2.3 » HYDEX 3.2.6 
-CDABS 3.2.4 »HYEXP 3.2.6 
-CDCOS 3.2.3 
-CDDIV 3.7 ~ I XXC 3.6 
- CDEXP 3.2.2 -IXXCD 3.6 
-CDLOG 3.2.4 -IXXI 3.6 
- CDMPY 3.7 ~IXXD 3.6 
- CDSIN 3.2.3 -IXXR 3.6 
-» CDSQRT 3.2.1 
- CDXXC 3.6 «RXXC 3.6 
- CDXXCD 3.6 -RXXCD 3.6 
- CDX XD 3.6 -RXXI 3.6 
-CDXXI 3.6 -RXXD 3.6 
- CDXXR 3.6 
- CEXP - 3.2.2 ~SIN 3.2.2, 3.2.3 
-CLOG 3.2.4 -SINH 3.2.3, 3.2.6 
-COS 3.2.2, 3.2.3 -SQRT 3.2.1, 3.2.5, 3.5.2 
~COSH 3.2.3, 3.2.6 
-CSIN 3.2.3 - TAN 3.2.3 
~CSORT 3.2.1 - TANH 3.2.6 
-CXDIV 3.7 
- CXMPY 3.7 
»CXXC 3.6 
~ CXXCD 3.6 
- CXXD 3.6 
- CXXR 3.6 
- CXXI 3.7 
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Function 


ABS 
ACOS 
A IMAG 
AINT 
ALOG 
ALOG1O 
AMAX 0 
AMAX1 
AMINO 
AMIN1L 
AMOD 
ANINT 
ASIN 
ATAN 
ATAN2 


CABS 
CCOS 
CDABS 
CDCOS 
CDEXP 
CDLOG 
CDSIN 
CDSORT 
CEXP 
CLOG 
C MP LX 
CONJG 
COS 
COSH 
CSIN 
CSQRT 


DABS 
DACOS 
DASIN 
DATAN 
DATAN2 
DBLE 
DCMPLX 
DCONJG 
DCOS 
DCO SH 
DDIM 
DEXP 
DFLOAT 
DIM 
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Math 
Function 


DIMAG 
DINT 
DLOG 
DLOG1O 
DMAX1 
DMIN1 
DMOD 
DNINT 
DPROD 
DREAL 
DSIGN 
DSIN 
DS INH 
DSQRT 
DTAN 
DTANH 


EXP 


FLOAT 


IABS 
IDIM 
IDINT 
IDNINT 
IFIX 
INT 
INT2 
ISIGN 


MAX0 
MAX1 
MINO 
MIN1 
MOD 


WNINT 
REAL 


SIGN 
SIN 

S INH 
SNGL 
SQRT 


TAN 
TANH 
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268 435 
536 870 
1 O73 741 
2 147 483 


4 294 967 
8 589 934 
17° (179) =S69 
34 359 738 


63 719 476 
137 488 953 
274 S77 906 
549 755 413 


1 099 511 627 
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ARITHNETIC 


TABLES 


TABLE OF POWERS OF TWO 
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TABLE OF POWERS OF SIXTEEN 


16 

4 

(Cts 

1 099 

17 592 

281 474 

4 505 599 

72 057 994 

1 152 921 504 


0 
1 
2 
3 
4 
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6 
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B 
Cc 
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HEXADECIMAL TO DECIMAL CONVERSION TABLE 


4,096 

3,192 
12, 288 
16, 384 
20,480 
24,576 
24,672 
32, 768 
36, 864 
40,960 
45,056 
49,152 
53, 248 
57,344 
61,440 


SMOADPOKENOUAS 


Decimal Values 


3, 072 
3, 328 
3, 584 
3,840 
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BYTE 
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TABLE OF MATHEMATICAL CONSTANTS 


Hexadecimal 


Decimal Value 
ecima ilue Value 


Constant 


49793 3.2431 


QSSNGL 43790 0.517C 


1.77245 38509 05516 1.C5BF 


1.14472 988548 49400 1.250D 


2. 71828 


18284 59045 2.B7E1 


0.36787 94411 71442 0.3E2D 


Ve 1.64872 12707 0012s 1. A612 
logy oe 0.43429 144819 03252 0.61 2D 
logge 1.44269 50408 448963 1.7154 
7 0.57721 56649 0.93C4 
Ln ¥ -0.54953 93129 $1645 -0.8CAE 
aa 1.41421 35623 73095 1. GA09 
Ln2 0.69314 71805 59945 0. B172 
logy 92 0.30102 99956- 63981 0.4D10 
V10— 3.16227 76601 68379 3. 29%B 
Ln10 2.30258 50929 94046 2.4D76 
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ASCLI CARD CODF CONVERS LON TABLE 


GRAPHIC GRAPHIC 


SPACE @ 40 8-4 
A 41 12-1 

8 B 42 12-2 

8 C 43 12-3 

8 D 44 12-4 

8 E 45 12-5 

2 F 46 12-6 

8-5 G 47 i2-% 

( 28 12-8-5 H 48 12-8 
) 29 11-8-5 I 49 12-9 
* 2A 11-8-4 J 4A 11-1 
+ 2B 12-8-6 K 4B 11-2 
; 26 0-8-3 L 4C 11-3 
= 2D 11 M 4D 11-4 
d 2E 12-83 N 4E 11-5 
/ 2F 0-1 O 4F 11-6 
0 30 0 P 50 11-7 
1 31 1 Q 51 11-8 
2 32 2 R 52 11-9 
3 33 3 S 53 0-2 
4 34 4 T 54 0-3 
5 35 5 U 55 0-4 
6 36 6 V 56 0-5 
7 37 q W 57 0-6 
8 38 8 Xx 58 0-7 
9 39 9 Y 59 0-8 
: 3A 8-2 Z 5A 0-9 
3B 11-8-6 [ 5B 12-8-2 
¢ 3C 12-8-4 = 5C 0-8-2 
= 3D 8-6 } 5D 12-8-7 
S 3E 0-8-6 ‘ 5E 11-8-7 
? 3F 0-8-7 se 5F 0-8-5 
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INDEX 


Absolute Error, 2-1, 3-2, 3-3, 37-5, 3-7, 3-9, 3-11 


Absolute Value, 2-1, 3-2, 3-5, 3-7, 3-9, 3-11, 3720, 3-21 


Absolute Value Function, 1-2, 3-19 
Accuracy of Function Subprograms, 2-1 
Accuracy, Related Terms, 2-1 

ANSI Standard X3J3, 1-1 

Argument, Complex, 3-3, 3-5, 3-7, 3-24 
Arithmetic, Floating-Point, 2-1 


Complex Argument, 3-3, 37-5, 3-7, 3-24 

Complex Arithmetic Function, 1-1, 172, 27-2, 3-24 
Complex Functions, 2-2, 3-4, 3-7 

Complex Number, 3-1, 3-3, 3-5, 3-19, 3-20, 3-24 
Complex Numbers, Manipulation of, 1-2, 3-19 
Complex Routines, 2~3/2-4 . 
Computation, Floating-Point, 2-1 

Conversions, Type, 27-2, 3716, 3-21 


Description of Function Subprograms, 3-1. 
Double~Precision Product of Single~Precision 
Factors Function, 1-2, 3-21 


Elementary Functions 
Exponential Function, 3-3, 3-4 
Hyperbolic Function, 3-11, 3-12, 3-13 
Inverse Trigonometric Function, 3-8, 3-9, 3-10, 3-11 
Logarithm Function, 3-7, 3-8 
Square Root Function, 3-1, 37-2 
Trigonometric Function, 3-4, 3-5, 3-6 
Error 
Absolute Error, 2-1, 3-2, 3-3, 3-5, 3-7, 3- 
Relative Error, 2-1, 3-2, 3-5, 3-7, 3-9, 3 
Error Messages, Generation of, 2-3/2-4 
Error Propagation, 2-1 
Exponential, 1-1, 2-2 
Exponential Function, 3-3, 3-4 
External Versions of Functions, 2-3/2-4 


Floating-Point Arithmetic, 2-1 
Floating-Point Computations, 2-1 
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Function 

Absolute Value Function, 1-2, 3-19 

Complex Arithmetic Function, 1-1, 3-24 

Double-Precision Product of Single-Precision 
Factors Function, 1-2, 3-21 

Exponential Function, 3-3, 3-4 

External Versions of Functions, 273/274 

Hyperbolic Functions, 3-11, 3712, 3-13 

Internal Versions of Functions, 2-3/2-4 

Inverse Trigonometric Functions, 3-8, 3-9, 3-10, 3-11 

Logarithm Function, 3-7, 3-8 

Manipulation of Complex Number Function, 1-2, 3-19 

Min and Max Functions, 1-2, 37-14, 3-15 

Positive Difference Function, 172, 37-20 

Remaindering Function, 1-2, 3-18 

Single Precision Function, 2-l 

Square Root Function, 3-1, 3-2 

Trigonometric Function, 3-4, 3-5, 3-6 

Truncation and Rounding Function, 1-2, 3-21 


Generation of Error Messages, 2-3/2-4 


Hyperbolic, 2-2 
Hyperbolic Function, 3-11, 3-12, 3-13 


Imaginary Part, 3-3, 3-5, 3-7, 3-24 

Internal Versions of Functions, 2-3/2-4 

Intrinsic Function, 2-3/2-4 

Intrinsic Version, 27-3/2-4, 3-1 

Inverse Trigonometric Functions, 3-8, 3-9, 3-10, 3-11 
Involution (Raising to Powers), 3-22 , 


Logarithm Function, 3-7, 3-8 


Manipulation of Complex Number Function, 1-2, 3-19 
Min and Max Functions, 1-2, 3-14 
Miscellaneous Functions 
Absolute Value Function, 1-2, 3-19 
Double-Precision Product of s:ingle~Precision 
Factors Function, 1-20, 3-21 
Manipulation of Complex Number Function, 1-2, 3-19 
Positive Difference Function, 1-2, 3-20 
Remaindering Function, 1-2, 3-18 
Transfer of Sign Function, 1-2, 3-20 
Truncation and Rounding Function, 1-2, 3-21 


Number, Complex, 3-1, 3-3, 3-5, 3-19, 3-24 


Overflow, 2-37 294; 3-4, 3-1:5,. 3718; 3-19, 3321; 3-22, 3-23, 3-24 
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Positive Difference Function, 1-2, 3-20 


Raising to Powers (Involution), 3-22 

Range Restrictions of Subprograms in the RTL, 2-2, 2-3 
Real Part, 3-3, 3-5 

Relative Accuracy, 2-1 

Relative Error, 2-1, 3-2, 3-3, 3-9, 3-ll, 3-21 
Remaindering Function, 1-2, 3-18 ; 
Routines, Complex, 2-3/2-4 


Single~Precision Functions, 2-1 
Square Root Function, 3-1, 3-2 
Subprograms 
Accuracy of Function Subprograms, 2-1 
Description of Function Subprograms, 3-1 
Range Restrictions of Subprograms in the RTL, 2-2, 2-3 
Transfer of Sign Function, 1-2, 3-20 
Type Conversion Subprogram, 1-2 


Terms Related to Accuracy, 2-1] 

Transfer of Sign Function, 1-2, 3-20 
Trigonometric Function, 3-4, 3-5, 3-6 . 
Truncation and Rounding Functions, 1-2, 3-21 
Type Conversion Subprogram, 1-2, 2-2 

Type Conversions, 2-2, 3-16, 3-21 


Underflow, 2-3/2-4, 3-23, 3-24 
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